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Amplified signal response in scale-free networks by collaborative signaling 
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Many natural and artificial two-states signaling devices are connected forming networks. The 
information-processing potential of these systems is usually related to the response to weak external 
signals. Here, using a network of overdamped bistable elements, we study the effect of a heteroge- 
neous complex topology on the signal response. The analysis of the problem in random scale-free 
networks, reveals that heterogeneity plays a crucial role in amplifying external signals. We have 
contrasted numerical simulations with analytical calculations in simplified topologies. 
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Signaling devices are the building blocks of many nat- 
ural and artificial information-processing systems, for ex- 
ample, cells in living organisms respond to their environ- 
ment by means of an interconnected network of receptors, 
messengers, protein kinases and other signaling molecules 
[l|, Q. In artificial systems, we remark recent studies of 
networks of fluxgate magnetometers, with potential ways 
to enhance the utility and sensitivity of a large class of 
nonlinear dynamic sensors e.g. the magnetometers, fer- 
roelectric detectors for electric fields, or piezoelectric de- 
tectors for acoustics applications by careful coupling and 
configuration [3|, |4|, l5| . Their proper functioning always 
imply high sensitivity to external signals. 

One of the most remarkable discoveries in non-linear 
physics, over the last thirty years, was the phenomenon 
of stochastic resonance. It represents the surprising ef- 
fect manifested in nonlinear dynamics where a weak sub- 
threshold input signal can be amplified by the assis- 
tance of noise. The requirements for such amplifica- 
tion are three: (i) a non-linear dynamical system endow- 
ing a potential with energetic activation barriers, (ii) a 
small amplitude (usually periodic) external signal, and 
(iii) a source of noise inherent or coupled to the system. 
Given these features, the response of the system under- 
goes resonance-like behavior as a function of the noise 
floor; hence the name stochastic resonance 0, 0j- 

The response of the system can be enhanced by cou- 
pling it through all-to-all or in arrays configurations 
[9( . This has been done mainly resorting to a linear cou- 
pling among oscillators. The most intriguing part of the 
phenomenon is that noise (a source of disorder) can fa- 
vor the amplification of a signal, an effect that usually 
requires a very precise synchrony (order). Recently in 
[llj . it has been shown that even in absence of noise a 
similar beneficial effects can be observed. There, different 
sources of diversity plays the role of the noise, inducing 
a resonant collective behavior. Moreover, it was pointed 
out that such a resonance effect does not depend on the 



regular network of connectivities. 

Our goal here is to give evidences that an equivalent 
amplification of a external signal can be obtained in a 
deterministic dynamical system of signaling devices in 
complex heterogeneous networks. 

The typical example, where the phenomenon has been 
investigated, consists on a bistable potential with a pe- 
riodic signal in a thermal bath. A general dimensionless 
equation describing this scenario is 



x - x 3 + Asin(ut) + rj(t) 



(1) 



where the bistable potential is V(x) — ~x 2 /2 + x 4 /i, 
the external signal Asm(ojt) and r?(i) a Gaussian white 
noise with zero mean and autocorrelation (rj(t)r)(t')) = 
2DS(t — t'). Eq.(p} represents the overdamped motion of 
a Brownian particle in a bistable potential and periodic 
forcing. In absence of forcing and noise, the system has 
two stable fixed points centered around ±1, which corre- 
sponds to the minimum of the potential energy function 
V{x). When the amplitude of the forced signal is sub- 
threshold, each node i oscillates around the minimum of 
its potential with the same frequency of the forcing sig- 
nal. 

The deterministic system we will study corresponds to 
a network of elements obeying Eq.([T]) in absence of noise. 
Noise has been explicitly excluded to avoid a possible un- 
controlled superposition of effects, that would make dif- 
ficult to determine the influence of the topology. The 
network is expressed by its adjacency matrix Aij with 
entries 1 if i is connected to j, and otherwise. For 
simplicity, from now on we consider only undirected un- 
weighted networks. Mathematically the system reads 



Asm(ujt) + XL 



i=l, 



• N, (2) 



where L; 



network, being k, 



is the Laplacian matrix of the 
the degree of node i. 



We have conducted numerical simulations of the sys- 
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Vi = (l-X)y l -yf + Asm{ujt) + Xx H , i = l,...,N (4) 
Now we make the following hypothesis: for a coupling A 
sufficiently small, the dynamics of the leaves can be de- 
coupled from that of the hub, obtaining then from Eq.Q 
the well known equation of an overdamped bistable oscil- 
lator. Moreover, when the signal is weak enough 
compared with the potential barrier, this can be 
conveniently linearized around one of the poten- 
tial minima. Then, Eq.([4]) can be solved for the ith 
node, and asymptotically for long time yields, 



FIG. 1: Average amplification (G) as a function of the cou- 
pling A, for both topologies, scale-free and all-to-all. Results 
corresponding to two different values of the average degree 
({k) = 3, marked with circles, and (k) = 5 with triangles) 
for the BA model network are plotted. Here N = 500, and 
u — 27rl0~ . The simulations were averaged over 1000 ran- 
domly chosen initial conditions. While the scale-free networks 
show a clear amplification of the external signal, for a signi- 
ficative range of values of A that depends on the average de- 
gree, the all-to-all connectivity does not amplifiy the signal. 



A 



[wcos(wi) — 2sin(wf)] (5) 



with £j = ±1 depending on the initial conditions. Insert- 
ing this solution into Eq.([3]) we obtain 



x H 



where 



-V H {x H ) + A 1 sm{ut) + B' cos(wt) + Xr]i (6) 



and all-to-all connectivity networks. While BA has be- 
come the paradigm of a growing model that pro vides with 
a scale-free (power-law) degree distribution [10(, all-to-all 
connectivity is a representative of homogeneous in degree 
networks. 

The results are depicted in Figure [T] Here the aver- 
age amplification (G) over different initial conditions has 
been computed as function of the coupling A, being the 
amplification defined as G = va.ax.iXi/ A. The frequency 
of the forcing signal u>, and the number of nodes N were 
here kept fixed. The obtained results show that in BA 
networks an amplification of the signal occurs, while in 
all-to-all networks no amplification of the external signal 
is observed for any value of the coupling A. 

This difference of behavior is attributed to the pres- 
ence of hubs in the BA networks, which are the main 
responsible for the heterogeneity in degree of the net- 
work. Keeping this in mind, we study a topology con- 
sisting on a star-like network (one hub and N-l periph- 
eral leaves). Such a topology is simple enough to be 
mathematically tractable, and simultaneously capable to 
capture the main trait of heterogeneity found in scale- 
free networks. Indeed, the resulting dynamical system 
decomposes in two parts: the dynamics of the hub (the 
highly connected node in the network with N — 1 links) , 
xh, and the leaves, j/j, linked to the hub. Then, the 
dynamical system @ becomes 



V H (x) = -[1 - X(N - l)]x 2 /2 + x A /A 



(7) 



x H = [1 - X(N - l)]x H 



r 3 

x H 



+ A sin(wt) + A yi (3) 



N-l 

£ 

i=l 



is the effective potential felt by the hub. The rest of 
parameters in Eq.([(j]) are r\i — X^i 1 ?^ anc ^ ^ = 
A[l + 2A(N-l)]/(w 2 +4), B' = -2X{N-l)/{u 2 +A). No- 
tice that, the problem has been reduced to the motion of 
an overdamped oscillator, in an effective potential driven 
by a reamplified forcing signal coming from the global 
sum of the leaves. The two possible solutions for the 
nodes Eq.©, namely oscillations around ±1, affect the 
dynamics of the hub as a quenched disorder represented 
by Xrji, and is equivalent to the diversity studied in 11]. 
Strictly speaking, the central limit theorem allows us to 
state that in the limit N — > oo, rji behaves as a random 
variable governed by a gaussian probability distribution 
with variance a 2 — N — 1 . This is so because the initial 
conditions were randomly chosen. The height of the po- 
tential barrier for the hub is now h = [1 — (N — l)A] 2 /4, 
and induces a strongly dependence on A via the factor 
(N — 1). Moreover, when A = 1/(N — 1) the barrier for 
the hub disappears, leading to a unique single fixed point 
Xh = 0. Note that both mechanisms, the quenched dis- 
order and the decrement of the potential barrier, rooted 
on the heterogeneity of the network, may cooperate to 
allow the hub to surmount the potential barrier. When 
A — 1), and for coupling and signal ampli- 

tude sufficiently small, Eq.© can be solved in the 
asymptotic long-time limit and yields 
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X H (t,T))t^oo ~ Xh ~ 



[{B'uj - A' a H ) sm(ui) - (B' a H + A'u) cos(wt)] 



'H 



(8) 
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FIG. 2: Average amplification (G) as a function of the cou- 
pling A. We compare the numerical simulation results of the 
dynamics of the star-like network Eqs, (13l4[l . and the theoret- 
ical result Eq. pop . Inset: Maximum average amplification 
{G)max as a function of the size of the network N. We com- 
pare again, the numerical simulation results of the dynamics 
of the star-like network Eqs. (|3l4 |l. and the theoretical approx- 
imation Eq, (|ll[) . Parameters are as in Fig. 1 



where an = (xjj ) being x^' the equilibrium points in 
absence of forcing. Note that xn(t,7]) depends also im- 
plicitly on the random variable n through the equilibrium 
point Xg to which the hub dynamics relaxes. Knowing 
the long time evolution for the hub, its amplification can 
be readily evaluated, and the result is 
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H 



A'a H ) 2 + (B' a H + A'uj) 2 
(9) 



In practice, we must average over initial conditions to 
cancel out the dependence on them. However, this turns 
out to be equivalent to averaging over 77, . Therefore, the 
average amplification is given by 



(G) 



1 



y/2ir(N - 1) J- 



drj e 



-G( V ) (10) 



The influence on the size of the network, N can be 
analyzed resorting to the equation above. In particu- 
lar, we are interested in the case of N large, and for 
A w 1/(N — 1), when the average amplification attains 
its maximum value. Making use of the Laplace method, 
it is possible to find an asymptotic solution to Eq.(JTUJ), 
in the limit N — > 00, which is 



(G) = G(0) + 1^3 



V + 0(N- 2 ), (11) 



FIG. 3: Maximum average amplification (G) m ax obtained for 
the star-like network and BA network as a function of the 
frequency of the external signal. The coupling is kept fixed 
to the value A = 0.04 for the BA network, and A = 0.002 for 
the star-like network. These coupling values are those that 
exhibit the maximum of (G) in Fig. 1 and Fig. 2, respectively, 
the solid line corresponds to the theoretical prediction for the 
star-like network given by Eo. (|10p . Other parameters are as 
in Fig. 1. 



Therefore, the maximum average amplification can be 
extracted from a network, turns out to be bounded from 
above by G(0). 

We have performed extensive numerical simulations for 
the dynamical system (O-Q, varying N and A. In Fig. [21 
the results of the amplification in Eq. (|10[) and the numer- 
ical simulations of the system equations are compared for 
different values of the coupling A, keeping fixed N and 
lu. The amplification grows monotonically with A until 
reaching a maximum value. Such a value coincides ap- 
proximately with A = l/(N — 1), which corresponds to 
the minimum value of the barrier height. Therefore, the 
maximum signal amplification occurs when the barrier 
of the effective potential of the hub attains its minimal 
value. This allows the hub to move back and forth be- 
tween the two minima of the effective potential, yielding 
a much wider oscillation. It is important to observe that 
a similar beneficial effect has been found in the stochastic 
resonance phenomena, but rather here it is due merely 
to the natural heterogeneity of the network. Then, it 
can be seen as a kind of resonance induced by topology, 
i.e. a sort of topological resonance. The inset shows the 
maximum amplification obtained for different values of 
N. Notice that for large N, the solution tends asymp- 
totically to a constant value as expected from Eq. lTPTj) . 
Finally, in Fig. [31 we compare the maximum amplifica- 



of the external signal. The amplification decreases mono- 
tonically with the frequency as predicted by the model. 

The remarkable agreement between simulations and 
analytical results validates our hypothesis of decoupling 
the hub from leaves. We are now in the situation of in- 
terpreting the results in the BA scale-free networks. A 
neat picture is given by considering that each highly con- 
nected node acts locally as it were the hub of a star-like 
network, with a degree k picked up from the degree dis- 
tribution. Therefore, for a given coupling A, we can find 
several star-like networks in different stages, depending 
on the degree of its local hubs. Recall that when A = 1/fc, 
the maximum signal amplification is attained for a hub 
with degree k. Furthermore, increasingly larger values of 
A would activate local hubs with smaller degrees. Since 
the network has several hubs, it should exist a wide range 
of values of A for which the maximum amplification is 
achieved following a cascade of amplification. In Fig. [T] 
such a behavior can be clearly observed. This contrasts 
with the results found for the star-like network, and an- 
ticipates that the scale-free network consists on a much 
more robust topology where such a new phenomena oc- 
curs. 

The mechanism described above remains valid until 
full synchronization occurs. For large values of A, the 
nodes become strongly coupled, and then they spatially 
synchronize behaving as a single node. We define the de- 
gree of spatial synchronization as the fraction of nodes 
\n + — n^\/N where n + and n_ are the number of ele- 
ments at the positive or negative well, respectively. In 
Fig. [5] we show how full synchronization emerges once a 
critical coupling is attained. In addition, numerical sim- 
ulations depicted in the same figure reveal that the path 
to synchronization is more pronounced for scale- free net- 
works with larger average degree. Then, for this case 
a smaller range of values for which maximum amplifica- 
tion is sustained should be observed in accordance with 
Fig. [T] Similar curves have been found in the prototype 
Kuramoto model [12j, which describes synchronization 
phenomena in nonlinear coupled oscillators, on top of 
complex networks [131 ]. However, the relevance of the 
present case stems from the fact that synchronization, 
contrary to which often happens in nonlinear coupled 
dynamics, has a negative effect in the amplification pro- 
cess. The all-to-all configuration in Fig. Q] presents full 
synchronization for any A different from 0. The reason 
for this observed behavior is that no source of disorder 
exist in this case. The spatial symmetry of the problem 
is unstable, and any small difference in concentration of 
elements in one of the two wells produces an asymmetry 
in the potential that evolves towards a monostable (syn- 
chronized) system. Note that this mechanism explains 
also the absence of amplification in Fig. [1] 

Summarizing, in this work we have shown that a topo- 
logical resonant-like effect emerges in scale-free topolo- 



Q- O BA seale-frce ( < k ) = 3 ) 
□■ ■□ BA scale-free ( ( k ) = 5 ) 

A A All-to-all 



,880000°°' 



,0° u 



X 



FIG. 4: Degree of synchronization of the BA and all-to-all 
networks versus the coupling. Two different average degree 
values for the BA network are plotted. The larger the average 
degree, the lower the critical coupling needed to achieve full 
synchronization is. The all-to-all network synchronizes for 
any value of A 7^ 0. Parameters are as in Fig. 1 



overdamped bistable dynamical system. The coopera- 
tive interaction of nodes connected to the hub due to 
the topology, has been shown to be the key feature for 
amplifying external signals. We have analytically proved 
that the phenomenon can be totally characterized in the 
framework of a simple topology consisting of a star-like 
network. The agreement between simulations and ana- 
lytical predictions allow us to understand the amplifica- 
tion curves as a function of the parameters of the system: 
the coupling, the size of the system and the frequency of 
the periodic external signal. At the light of the current 
results, we speculate that the ubiquity of scale-free topo- 
logical structures in biological systems can be related to 
its ability to amplify sensitivity to weak external signals. 
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